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Abstract 



We review diffractive deep inelastic scattering (DIS) in the light of the coUincar factor- 
ization theorem. This theorem allows to define diffractive parton distributions in the leading 
twist approach. Due to its selective final states, diffractive DIS offers interesting insight into 
the form of the diffractive parton distributions which we explore with the help of the satura- 
tion model. We find Regge-like factorization with the correct energy dependence measured at 
HERA. A remarkable feature of diffractive DIS is the dominance of the twist-4 contribution 
for small diffractive masses. We quantify this effect and make a comparison with the data. 



1 Introduction 



A significant fraction (around 10%) of deep inelastic scattering (DIS) events observed at HERA 
at small-x are diffractive events The proton in these events escapes almost unscattered 

down the beam pipe, losing only a small fraction xjp of its initial momentum. The slightly 
scattered proton is well separated in rapidity from the rest of the scattered system, forming a 
large rapidity gap, the characteristic feature of diffractive DIS. The ratio of diffractive to all DIS 
events is to a good approximation constant as a function of Bjorken-j; and Q^. Thus in a first 
approximation, DIS diffraction is a leading twist effect with logarithmic scaling violation. For 
recent reviews on diffractive DIS see |5| . 

Historically, the first description of diffractive DIS was provided in terms of the Ingelman- 
Schlein (IS) model |^]. The model is based on Regge theory in which diffractive processes are due 
to the exchange of a soft pomeron. The novelty of the IS model lies in the assumption that the 
pomeron has a partonic structure as do real hadrons. The diffractive structure function factorizes 
into a "pomeron flux" and a pomeron structure function. The latter function is written in terms 
of the pomeron parton distributions which obey the standard DGLAP evolution equations. The 
initial conditions for these equations are determined from a fit to the data ||l], |2|, or using the 
phenomenology of soft hadronic reactions [p. Despite conceptual difficulties (the pomeron is 
not a particle) this idea turned out to be very fruitful in the description of the data. 

The IS model brings the issue of coUinear factorization into the leading twist description 
of DIS diffraction. By this we mean the consistent factorization of the diffractive cross sec- 
tion into a convolution of diffractive parton distributions which satisfy the DGLAP evolution 
equations and computed in perturbative QCD hard cross sections, in full analogy to inclusive 
DIS P, |l^]. Collinear factorization has been rigorously proven by Collins for diffractive 
DIS . Factorization, however, fails for hard processes in diffractive hadron-hadron scattering 



|13|. The IS model assumes collinear factorization, imposing an additional assumption on the 
xjp-dependence of the diffractive parton distributions (called Regge factorization). In the gen- 
eral framework of collinear factorization the diffractive final state is treated fully inclusive, in 
particular, the mechanism leading to diffraction is not elucidated. The IS model is an attempt 
to provide such a mechanism. 

The detailed description of diffractive processes in DIS, starting from perturbative QCD, is 
achieved by modelling the diffractive final state as well as the interaction with the proton. Such 
an analysis goes beyond the leading twist description. In a first approximation, the diffractive 
system (separated by the rapidity gap from the proton) is formed by a quark-antiquark [qq) 
pair in the color singlet state A higher order contribution is represented by the qq 

pair with an additional gluon g emitted [16, 17, In the simplest case, the colorless 

exchange responsible for the rapidity gap is modelled by the exchange of two gluons coupled 



to the proton with some formfactor [20, 21 [ or to a heavy onium which serves as a model of 
the proton |2^. Higher order corrections are included by the the BFKL summation of gluon 
ladders [|2^ or using the color dipole approach [p3. The diffractive processes in DIS has also 



been studied within the semiclassical approach in |25|[. In particular, the relation between the 



physical pictures in different frames was elucidated in [26[. Let us add that a simple physical 
interpretation of diffractive scattering emerges in the proton rest frame, where the formation 
of the diffractive system is stretched in time. Far upstream the beam pipe the virtual photon 
decays into a virtual diffractive system which then elastically scatters on the proton (without 
the color exchange), picking up energy to form a real state. 
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The immediate problem faced in such a modelhng is the strong sensitivity to nonperturbative 
effects due to the dominance of aligned jet configuration (discussed in Section |^. Thus, we need 
a description of the interactions in the soft regime. The saturation model [^^, which has 
already been very successful in describing inclusive and diffractive DIS data [E^], provides such 



a description (see also ^ ^ for related analyses). Recent theoretical studies Q 
justify the assumed analytic form of this model. 

In this paper we are going to address questions related to collinear factorization in the above 
approach to DIS diffraction. How to find the diffractive parton distributions? What is their 
form? Do they support Regge factorization? Is the leading twist contribution sufficient in the 
description of DIS diffraction? How important are higher twist contributions? Answers to all 
these questions will be found by modeling the diffractive system with the help of perturbative 
QCD. The interaction with the proton is described by the saturation model |27|. From this 
perspective, we critically reexamine the assumptions made in the Ingelman-Schlein model. 

Summarizing our results, we find the diffractive quark and gluon distributions which serve 
as the initial conditions for the DGLAP evolution equations. Due to a specific form of the satu- 
ration model, the Regge type factorization in xjp is found, although Regge theory has not been 
applied. Moreover, the correct energy dependence of diffractive DIS measured at HERA is ob- 
tained. We also perform the numerical analysis for the comparison with the diffractive data. As 
expected, see iQ, the twist-4 contribution from the gg-pair produced by longitudinally polarized 
photons plays a crucial role in the region of small diffractive mass (M <C Q). The leading twist 
description with the DGLAP evolution is insufficient in this regime and the twist-4 component 
(suppressed by 1/Q^) accounts for the difference between the leading twist contribution and the 
data. Thus, there is no need in our analysis for a singular (or strongly concentrated at /3 ~ 1) 
gluon distribution as in the pure leading twist description Q. The universality of the energy 
dependence, assumed in the IS model, is broken by the twist-4 contribution which has a steeper 
xjp-dependence than the leading twist contribution. The first indication of this effect seems to 
be observed in the data 0|. 

The paper is organized as follows. In Section |2| we give a general introduction to diffractive 
parton distributions. In Section ^ we present the initial distribution which we extracted from 
the perturbative QCD approach. In Section |^ we make a comment on Regge factorization in 
the context of the saturation model. Numerical results are presented in Section ^ and we finally 
conclude with a brief summary in Section |^. 



2 Diffractive parton distributions 

There are two dimensionful variables, the mass of the diffractive system M^, and the momentum 
transfer t = {p — p')'^, which characterize diffractive DIS (see Fig. |^). They come in addition 
to and W'^ which are well known from inclusive DIS. 

The mass and energy related variables, and W'^, are usually rewritten in terms of the 
following dimensionless variables: 

Q^ + M^-t 

which describes the fraction of the incident momentum lost by the proton or carried by the 
pomeron, and 
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Figure 1: Kinematic invariants in DIS diffraction. 



which is the Bjorken variable normahzed to the pomeron momentum. The true Bjorken variable 
X connects the two variables ^ 

^ = Q2TW^ = 

In the following we neglect t in the definition of xjp and /? since usually |t| ^ Q^, M^. 

The diffractive structure functions depend on four invariant variables (x, Q'^,X]p, t), and their 
definition follows the inclusive counterparts by relating to the diffractive DIS cross section 



l + n-yf] AIL - y2 AIL] , (4) 
dx dQ"^ dxjpdt x 1 L i dxjpdt dxpdt j 

For simplicity, we introduce the following notation 

F2^i^)(x,Q2,xj.,t) ^ ^lL^x^Q2^xjp^t). (5) 

As usual, -F2^^^'* = F^^^^ +F^^^\ where T and L refer to the polarization of the virtual photon, 
transverse and longitudinal, respectively. The structure function F^^^^ has the dimension GeV~'^ 
because of the differential dt in the cross section (^). The dimensionless structure function F^(^) 
is defined by integrating F^^^) over t, 

F^}l\x,Q^xip) = J dtFp^\x,Q^,xip,t). (6) 

^D(3) jg nieasured when the final state proton transverse momentum is not detected. 

The diffractive parton distributions ^^/p introduced according to the collinear factoriza- 
tion formula Q, 

F^^^\x,Q\xjp,t) = n d^T^i^{i,^?,xjp,t)Ca{x/i,Q^/f?), (7) 

with a = q,g denoting a quark or gluon distribution in the proton, respectively. In the infinite 
momentum frame the diffractive parton distributions describe the probability to find a parton 
with the fraction ^ of the proton momentum, provided the proton stays intact and loses only 
a small fraction xjp of its original momentum. Ca are the coefficient functions describing hard 
scattering of the virtual photon on a parton a. They are identical to the coefficient functions 
known from inclusive DIS, 

Caix/t QV/^') = el 6{1 - x/0 + Oias) . (8) 
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Formula (0) is the analogue of the inclusive leading twist description for inclusive DIS. 
The inclusive structure function F2 is factorized in a similar way into computed in pQCD 
coefficient functions and nonperturbative parton distributions. The scale //^ is the factor- 
ization/renormalization scale. Since the l.h.s of eq. (|^) does not depend on this scale, i.e. 

we find the renormalization group equations (evolution equations) for the 
diffractive parton distribution 

Ai^^ •^i/p(C,/^^a^JP,^) = — Pa/b{^/z,as{lJ^)) J'kip{z,^?,X]p,t) , (9) 

^ b ^ 

where Pa/b are the standard Altarelli-Parisi splitting functions in leading or next-to-leading 
logarithmic approximation. Since the scale fj, is arbitrary, we can choose fi = Q ^ Aqcd- With 
this scale the evolution equations are usually presented. 

The integration in and @ is only done up to the fraction xjp of the proton momentum, 
since the active parton cannot carry more than this fraction of momentum. The proton remnants 
carry the remaining fraction (1 — xjp). If we refer the longitudinal momenta of the partons to 
X]pp instead of the proton total momentum p, the structure functions and parton distributions 
become functions of j3 = x/xjp or /?' = i/xp. With this notation, we rewrite (0) and (P) in the 
following form: 

F^^^\p,Q\xjp,t) = J2 I'dp' xjpT^/piP',ti^xjp,t) CaiP/f3',Q^/ti^) (10) 



and 



d dz 
f^^^J'^/p{P,f^\xjp,t) =11]^ —Pa/b{P/z,as{fi')) Ti;;jX^,fi^xiP,t). (11) 

Thus, we obtain a description similar to inclusive DIS but modified by the additional variables 
X]p and t. Moreover, the Bjorken variable x is replaced by its diffractive analogue (3, eq. (§). 
Notice that x jp and t play the role of parameters of the evolution equations and does not affect 
the evolution. According to the factorization theorem the evolution equations (|l^) are applicable 
to all orders in perturbation theory. 

In the lowest order approximation for the coefficient functions (^), we find for the diffractive 
structure function 

F^^^\l3,Q\xjp,t) = J2 elPxjpT^/piP,Q^xjp,t), (12) 

a=q,q 

where the sum over the quark flavours is performed. 

The collinear factorization formula ( p!0| ) holds to all orders in for diffractive DIS |12]. 
However, this is no longer true in hadron-hadron hard diffractive scattering [^, ^], where 
collinear factorization fails due to flnal state soft interactions. Thus, unlike inclusive scattering, 
the diffractive parton distributions are no universal quantities. The can safely be used, however, 
to describe hard diffractive processes involving leptons. A systematic approach to diffractive 



parton distributions, based on quark and gluon operators, is given in [11, |2^ . 

Until now, our discussion has been quite general, in particular we have not referred to the 
pomeron. In the Ingelman-Schlein (IS) model diffraction is described with the help of the 
concept of the soft pomeron exchange. In addition, it is assumed that the pomeron has a hard 
structure. In DIS diffraction, this structure is resolved by the virtual photon, as in the standard 
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DIS processes. Following the results of Regge theory, the IS model is based on the assumption 
of Regge factorization. In the context of the diffractive parton distributions it means that the 
following factorization holds 11| 



XjpJ^^/j,{P,Q^,Xjp,t) = f{xjp,t) fa/AP,Q^), (13) 
where the "pomeron flux" f{x]p,t) is given by 

/(x^,O = 0xi.--W. (14) 

Thus, the variables {x]p,t), related to the loosely scattered proton, are factorized from the 
variables characterizing the diffractive system Q^). B{t) is the Dirac electromagnetic form 
factor ||3^ , ajp(t) = 1.1 + 0.25 GeV~^-t is the soft pomeron trajectory |^] and the normalization 
of f{x]p,t) follows the convention of |Q. The function fa/]p{P,Q'^) in eq. ( [l^ ) describes the 
hard structure in DIS diffraction, and is interpreted as the pomeron parton distribution. Now, 
the diffractive structure function (|l^) becomes 

F^^^\p,Q^xjp,t) = fixjp,t) elPfa/AP,Q^), (15) 

a=q,q 

where the summation over quarks and antiquarks is performed. The Q^-evolution of fa/ipiP, Q"^) 
is given by the DGLAP equations ([ll|). The t-dependence in the pomeron parton distributions 
is neglected. 

The pomeron parton distributions are determined as the parton distributions of real hadrons. 
Some functional form with several parameters is assumed at an initial scale and then the param- 
eters are found from a fit to data ||l], ^, ^ using the DGLAP evolution equations. An alternative 
approach for the determination of the initial distributions makes use of the phenomenology of 
soft hadronic reactions 

Ingelman and Schlein have conceived of their approach primarily to describe diffractive hard 
scattering at hadron colliders. This includes the concept of Regge factorization as well as 
pomeron parton distributions. Unfortunately, collinear factorization was proven to be wrong in 
the case of hard diffraction at hadron colliders. Although it holds for diffractive DIS, this does 
not imply Regge factorization nor the existence of pomeron parton distributions. Our approach 
discussed below does not make use of either Regge factorization or pomeron parton distributions. 
It does, however, result in diffractive parton distributions and xjp-factorization. 



3 Diffractive parton distributions and the saturation model 

In ref. an analytic form for the dipole cross section a{x, r) was suggested, based on the idea 



of saturation |35], which allows to describe the proton structure function F2 at small x, 

F2{x,Q^) = -^- I Svdz (\^T{r,z,Q^)\^ + \^L{r,z,Q^)\^) a{x,r), (16) 

47r Oerra J ^ ' 

where \I't,l is the qq dipole wave function for transverse (T) and longitudinal (L) photons, r 
is the dipole transverse size and z is a fraction of the virtual photon momentum carried by a 
quark (antiquark). A distinctive feature of the dipole cross section is its scaling form, i.e. 

a(x,r) = g{r I Rfi(x)) = (Jq |i - exp(-r^/i?o(x)| , (17) 



5 



Figure 2: The diffractive qq and qqg contributions to F,^^^\ 

where the function Ro{x) = 1/GeV (x/xq)^/'^ monotonically vanishes when x ^ 0, and ctq is an 
overah normahzation. The three parameters of erg = 23 mb, xq = 3 ■ 10"^ and A = 0.29 were 
found form a fit to inclusive DIS data at small-x. 

The postulated form allows the transition to saturation where 6" = cq for large dipole sizes, 
while for small r color transparency is assumed, cj ~ r^. The fact that the dipole cross section 
depends on r and x through the dimensionless ratio r/Ro{x) leads to the prediction of a new 
scaling property for small-x data ||39t| . 

The presented model successfully describes F2 at small x, including the region of small 
values. With regard to diffractive processes in DIS, a leads to the same dependence on x and 

of diffractive cross section as for inclusive DIS, and gives a good description of the data 
p8|| . Thus, the constant ratio between the diffractive and inclusive cross sections finds a natural 
explanation in this model. 

Following the idea of the analysis the diffractive structure function F^^^^ is the sum 
of the three contributions shown in Fig. |2|, the qq production from transverse and longitudinal 
photons, and the qqg production, 

F^^''\(3,Q^xjp) = F];^ + F,^- + f4 , (18) 

where T and L refer to the polarization of the virtual photon. For the qqg contribution only the 
transverse polarization is considered, since the longitudinal counterpart has no leading logarithm 
in . In this approach, the diffractive qq and qqg systems interact with the proton like in the 
two gluon exchange model. The coupling of the two gluons to the proton is effectively described 
the dipole cross section, determined from the analysis of inclusive DIS. 



The computation of the diffractive structure functions in this case was presented in [28]. 
Here we quote only the final results. The transverse qq part is given by 

-, 2/? A:2 



™ Mtt'^Bdxp ^ ■> [l-py Jo 1^ 4/3 A;2 
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and the longitudinal qq contribution takes the form 



where the "impact factors" (/)j=o i read 




(^i(fc,/3,xjp) = li' / <irrK,( ^/^^fcrj Ji(A:r) <7(xjp, r) , (21) 

and Ki and Jj are the Bessel functions. The quoted formulae correspond to eqs. (32) and (33) 
in I^S), respectively, with the angular integration already done, see also the second paper of 
p9| for a similar result. is the diffractive slope, present due to the integration over t of 
F^(^) with the assumed t-dependence: F^^^^ = exp{B£)t) F^^^\ Its value is taken from the 
HERA experiments. The A;- integration in eqs. ( [l9| ) and (^0|) is performed over the transverse 
momentum of the quarks in the qq pair, k = |k_L|. The final result depends on the squared 
"impact factors", and thus on the square of the dipole cross section a. 

One should realize that the longitudinal contribution is suppressed by a power of in 
comparison to the transverse contribution, i.e. the longitudinal structure function is a higher 
twist contribution. Although of higher twist nature this contribution has some importance as 
will be seen later. 

The leading twist part of the transverse structure function, which corresponds to the diffrac- 
tive qq production, can be extracted from (|^) by neglecting the factors with powers of k'^/Q'^ 
under the integral and taking the upper limit of the integration to infinity. Strictly speaking, 
energy conservation is violated in this case, but the corrections are of higher twist nature which 
in this case will be neglected because of their smallness. With the new limit the integral is still 
finite, and we can write the leading twist part of (|Tg|) as 

- uii^ Y'lSwC ^■<^- ■ <^^' 

A detailed analysis based on the dipole representation of the qq wave function shows that the 
approximation leading to eq. ( p2| ) corresponds to the aligned jet configuration of the qq pair in 
the proton rest frame. The smallness of the factor 



1 - /? g2 j\^2 

which we neglected in eq. ( [l9| ) [M is the diffractive mass and z, (1 — z) are the longitudinal 
momentum fractions of the final state quarks with respect to the photon momentum) means 
that one of the quarks takes almost the whole longitudinal photon momentum (e.g. z ~ 1) while 
the other quark forms the remnant with {\ — z) ~ 0. Similar conclusions on the leading twist 



part of diffractive DIS have been drawn in ref. |26|. 



Now, we can determine the diffractive quark distributions according to eq. (|l2|) (integrated 
over t), 

^5?^^^ = 2 E 4 l""^^^ Q'^^^JP) (24) 
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Hence the dijjractive quark distribution is given by (independent of the quark flavour /) 



dk^ cj)l{k,(3,xip), (25) 



Notice the lack of the Q^-dependence on the r.h.s. of eq. (|2^). This may be viewed as a 
consequence of not having included ultraviolet divergent corrections which would require a cutoff. 
With those corrections the parton distributions become ^^-dependent, and evolution would relate 
the distributions at different ^'^ values. Still, we may use the found diffractive quark distributions 
as the input distributions for the DGLAP evolution equations at some initial scale. Of course, 
the choice of the initial scale introduces an uncertainty for the prediction. 

The detailed discussion of the qqg contribution can be found in (2^ and with the details 
of the calculations in the Appendix. The new contribution was computed assuming strong 
ordering in transverse momenta of the gluon and the qq pair, i.e. k±g <C k±q ~ k±q. This 
assumption allows to treat the qqg system as a gg dipole in the transverse configuration space 
r, where r is the Fourier conjugate variable to the quark transverse momentum k^. 

The formula for the qqg diffractive structure function which we quote below corresponds to 
eq. (39) in [28|, integrated over the azimuth angle in configuration space Q. Thus we have 



F^ggiP,Q\xjp) 



81/3 



2 ^ 

256tt^BdX]p 27r Jp (3' 



dl3' 



1 



+ 



^x2- 



(26) 



{l-P'f Jo 



where the impact factor 02 is given by eq. ( ^l|) with i = 2. The variable (i' describes the 
momentum fraction of the i-channel exchanged gluon with respect to the pomeron momentum 
xjpp. The combination k'^ /{I — f3') which enters the logarithm is its mean virtuality, and k = |k_L| 
is the transverse momentum of the final state gluon. The term in square brackets under the first 
integral is the Altarelli-Parisi splitting function for g — s- qq, which results from the approximation 
that the transverse momentum of the emitted gluon is smaller than transverse momenta of the 
quarks. 

The diffractive gluon distribution can be found from eq. (^6|). In the calculation of this 
contribution strong ordering between the gluon and quark transverse momenta was assumed. 
In this approximation the integration over the transverse momentum of the quark loop gives a 
logarithmic contribution which has a natural lower cutoff, the virtuality of the gluon k"^ / {1 — j3'). 
At the same time the virtuality of the gluon should not exceed Q^. This is the origin of the 
logarithmic term in (^6|). Collinear factorization means that we can pull that logarithm out of 
the integral over the gluon transverse momenta, and add to it an arbitrary initial scale Qo < Q^- 
Thus we can write 



QQg 



2E 

/ 



a. 



2tt 



1 (0_ 

■ T 



/3V V/3' 



X]pJ^g/p{P',X]p) 



(27) 



where the diffractive gluon distribution is given by 



g^(j3,X]p) = X]pJ^^/p{f3,X]p 



■D 



81 



2bQTT^BDXip (1 - I3Y 



dk"^ (pl{k,(3,xjp) (28) 



^ A factor 1/2 missing in eq.(39) of ref. |28| was correctly pointed out in ref. 
numerical results in [ESI. 



This does not affect the 
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and (j)2 is given by eq. ( pl| ) with i = 2. As in the case of the quark distribution (p5|), the found 
gluon distribution does not depend on Q^, and serves as the initial distributions at some fixed 
scale Qq. 

The motivation for the above identification of the diffractive gluon distributions is the struc- 
ture in the curly brackets on the r.h.s of eq. (|27|). It is identical to the structure resulting from 
the DGLAP evolution with one splitting of the gluon into a qq pair. 

The combined initial parton distributions (^) and ( [2^ ) ( depicted in Fig.^ ) allow a complete 
description of the leading twist part of diffractive DIS by serving as the initial conditions for 
the DGLAP evolution equations. DGLAP evolution means that our diffractive system becomes 
more complicated due to additional parton emissions. 

The longitudinal, higher twist contribution requires a separate treatment. It becomes im- 
portant for large values of P, where the qq and the qqg production from transverse photons is 



negligible |28, 16, 15]. In our present analysis we simply add this contribution to the evolved 



leading twist part. The complete expression of the structure function reads 



where F^^^^^'^ ig given by 



pm) ^ F^^{3)(^2^) + (29) 



^D(3)(LT) ^ 2^:6^/3 (A g^:^^), (30) 
/ 



with the full DGLAP evolution. F^q is given by eq. 

In the following we present a comparison of the description based on eq. ( p^ ) with the 
diffractive data from HERA. Before doing that we briefly discuss xjp-factorization. 



4 The issue of Regge factorization 

Up to now, we have not made use of the particular form (|l7|) of the dipole cross section a. We 
only implicitly assumed that the integrals involving a are finite. The scaling property, i.e. that 
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fj is a function of the dimensionless ratio r/RQ(x), has the remarkable consequence that the 
xjp-dependent part of the found diffractive parton distributions (DPD) factorizes. 

Introducing the dimensionless variables k = kRQ{x) and f = r/RQ(x) in (^) and (^) and 
assuming Qq to be a fixed scale, we find the following factorized form 

We have introduced a notation similar to that in ([T^ for the /3-dependent factors. This type of 
factorization looks like Regge factorization but has nothing to do with Regge theory. It merely 
results from the scaling properties of the saturating cross section a. Since the evolution does 
not affect the xjp-dependence of the DPD, the factorized form will be valid for any scale Q^. 

Now, we can rewrite eq. pO|) as 



(33) 

in which the xjp-dependence is factored out. In the saturation model {Ro{x) ~ x^/^) the 



parameter A = 0.29 was determined from a fit to inclusive DIS data only |27]. The same value 
holds for diffractive interactions, thus we find a definite prediction for the xjp-dependence of the 
leading twist diffractive structure function 

^ ^^i-A . (34) 



At present, the bulk of diffractive data in DIS support the factorized form (plj). They are 
usually interpreted || in terms of the t-averaged pomeron intercept cFp, i.e. 

^2^(3) ^ ^^-2^^ (35) 

Such a dependence has been introduced in the spirit of the Ingelman-Schlein model (|l5|), with 
the t-integration performed, F^^^^ ~ / dt f{xip,t) ~ x^^"^. Thus, according to and (|35|) 
we find 

OF = ^ + 1 « 1.15, (36) 

which is in remarkable agreement with the values found at HERA, cFp = ajp(O) — 0.03 = 1.17 
by HI g and oF = 1.13 by ZEUS @. 

Summarizing, the leading twist description extracted from the saturation model of DIS 
diffraction leads to the factorization of the xjp-dependent part of the cross section similar to 
Regge factorization. It correctly predicts the value of the "effective pomeron intercept" . The 
Q'^-dependence of the diffractive structure function does not affect the xjp-factorization. This 
means that the saturation model for the dipole cross section gives effectively the result which 
coincides with the Regge approach to DIS diffraction, although the physics behind is completely 
different. The relative hardness of the intrinsic scale 1/Rq{xip) ~ 1 GeV in the saturation model 
suggests that DIS diffraction is a semihard process rather than a soft process as Regge theory 
would require. 
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In the presented description, the leading twist structure function vanishes when /3 — > 1, i.e. 
for small diffractive mass case for the higher twist longitudinal 

contribution Fj^^^, eq. (20), which dominates in the region of /3 ~ 1 The expected xp- 



dependence for this contribution is given by 

X]pRl{x]p) 

which clearly violates the universality of the effective pomeron intercept, assumed in the Ingelman- 
Schlein model. The first indication of this effect is observed at HERA by measuring the effective 
pomeron intercept in different regions of the diffractive mass as a function of Q. The in- 
tercept seems to be bigger for smaller diffractive masses (/3 1), and the description based on 
perturbative QCD gives a satisfactory explanation 



5 Comparison with data 

In this section we present a comparison of the leading twist part of diffractive DIS including 
DGLAP evolution plus the corresponding longitudinal twist-4 component (eq. (p9|)) with the 
diffractive data from the HERA experiments. 

In Fig. ^ we show the quark and gluon diffractive distributions at some initial scale Qo- This 
scale, however, is not determined in this approach. Thus, it can be treated as a phenomenological 
parameter which may be tuned to obtain the best description of the DIS diffractive data. We 
tried various options and found that Qg ~ l-^ GeV^ is the best choice for this purpose^. We 
used the leading logarithmic evolution equations with three massless flavours, and the value of 
Kqcd = 200 MeV in a,. 

Fig. 1^ shows the results of our studies with data from HI and Fig. ^ data from ZEUS. 
The dashed lines represent the total contribution according to eq. (pQ]). The pure leading twist 
structure function F2 with the leading DGLAP evolution is shown by the dotted lines. 

The difference (if visible) between the dashed and dotted lines is the effect of the longitudinal 
twist-4 component F^^ added into the leading twist result. As expected, the twist-4 component 
is significant in the large-/? domain, see also for a detailed discussion. Notice the change 
of the slope in xjp when the twist-4 component is added. The overall agreement between the 
data and the model ( [2^ ) is reasonably good, taking into account the fact that the only tuned 
parameter is the initial scale for the evolution Qo- The parameters of the evolution equations 
are standard, and the diffractive slope Bd = 6 GeV~^ which effects the normalization is taken 
from the measurements as in the analysis [^8|. The significance of the twist-4 contribution at 
large (3 is clearly demonstrated. 

Looking at Fig. |^ and ^, we realize that at large /3 and the agreement between the data 
and the model with evolution starts to deteriorate. The reason for this is illustrated in Fig. ^ 
where we show F2' for fixed value of rcjp = 0.0042 as a function of /3, for different values of 
Q^. The DGLAP evolution depopulates the region of large f3, shifting the parton distributions 
towards smaller values of (3 (compare the dotted lines showing the initial distributions and the 
dashed lines showing the evolved distributions). The twist-4 component F^^ (the difference 
between the solid and the dashed lines) largely compensates for this effect. Its significance, 
however, diminishes as rises due to the l/Q^ dependence of twist-4. 

■^This value has been changed in comparison to the journal version because of a numerical error in the analysis 
there. The results of a comparison with the data are practically the same. 
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There are at least two effects which may enhance F^^ when rises, thus accounting for 
the difference between the data and the discussed description. In principle, twist-4 should also 
have logarithmic evolution, in addition to the 1/Q^ dependence, which could push F^^ into the 
right direction. Unfortunately, this aspect is beyond the scope of the present paper. Another 
effect which is important at large /3 is skewedness of the diffractive parton distributions, see |40|] 

T L — 

for a recent discussion and references therein. This effect is known to enhance F^q at large /3. 
Since the enhancement rises with Q^, the skewedness may account for the discussed difference 
between the data and their description in the region of large (3 and Q^. 



6 Conclusions 



We have reviewed the description of diffractive deep inelastic scattering in the light of the 
collinear factorization theorem. This theorem applies to the leading twist terms of the cross 
section and introduces the notion of diffractive parton distributions. We have extracted a precise 
analytic form for these distributions from the approach in which the diffractive state is formed 
by the qq and qqg systems, computed in pQCD. The convolution with a dipole cross section from 
the saturation model leads to the xjp-factorization of the diffractive structure functions, similar 
to Regge factorization, which correctly describe the energy dependence found at HERA by the 
HI and ZEUS collaborations. We further evolved the diffractive parton distributions with the 
DGLAP evolution equations and pointed out the significance of the twist-4 component at large 
(3 for the agreement with the data. The latter was originally advocated in [l^ and originates 
from the longitudinal qq contribution. The twist-4 component breaks the universality of the 
effective xjp-dependence, making it stronger at large (3. This stays in contrast to the assumed 
universality in the Ingelman-Schlein model. Finally, we suggest some necessary modifications 
of the description at large (5 and to improve the agreement with the data in this kinematic 
region. 
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H1 DATA (1994) 
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Figure 4: The comparison with HI data The dashed hnes correspond to the leading twist 
contribution with the twist-4 component added, eq. (|29|) . The leading twist contribution is 
shown by the dotted lines. 
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ZEUS DATA (1994) 
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Figure 5: The comparison with ZEUS data Q. The dashed hnes correspond to the leading 
twist contribution with the twist-4 component added, eq. (p9|). The leading twist contribution 
is shown by the dotted lines. 
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Figure 6: The effect of the DGLAP evolution in . The dotted hues show the leading twist 
structure function at the initial scale Q\ = "i GeV^. The dashed lines correspond to evolved 
structure function at the indicated values of . The solid line is the sum of the evolved leading 
twist contribution and twist-4 component, eq. p9|). The ZEUS data at xjp = 0.0042 are shown. 
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